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ABSTRACT 

We develop a new "Multicomponent and Variable Velocity" (MVV) galactic outflow 
model for cosmological smoothed particle hydrodynamic (SPH) simulations. The MVV 
wind model reflects the fact that the wind material can arise from different phases in 
the interstellar medium (ISM), and the mass- loading factor in the MVV model is a 
function of galaxy stellar mass. We flnd that the simulation with the MVV outflow has 
the following characteristics: (i) the intergalactic medium (IGM) is hardly heated up, 
and the mean IGM temperature is almost the same as in the no-wind run; (ii) it has 
lower cosmic star formation rates (SFRs) compared to the no-wind run, but higher 
SFRs than the constant velocity wind run; (iii) it roughly agrees with the observed 
IGM metallicity, and roughly follows the observed evolution of Vl(Qj iv); (iv) the lower 
mass galaxies have larger mass- loading factors, and the low-mass end of galaxy stellar 
mass function is flatter than in the previous simulations. Therefore, the MVV outflow 
model mildly alleviates the problem of too steep galaxy stellar mass function seen 
in the previous SPH simulations. In summary, the new MVV outflow model shows 
reasonable agreement with observations, and gives better results than the constant 
velocity wind model. 

Key words: method : numerical — galaxies : evolution — galaxies : formation — 
galaxies : high redshift — galaxies : mass function — cosmology : theory 



1 INTRODUCTION 

It has been pointed out by many researchers that galac- 
tic outflow plays a crucial role in galaxy formation. For ex- 
ample, the physical size and mass of a disk galaxy might 
be deter mined by the bala n ce between gas accre tion and 
outflow (iDekel Silklll986l : Iwhite h FrenkI Il99lh . Galac- 
tic outflows also regulate the star formation by removing 
gas from galaxies, and help to alleviate the "overcooling" 
problem seen in the earlier simulations of galaxy formation. 
In addition, the observations of quasar absorption lines re- 
veal that the intergalactic medium (IGM) is considerably 
enriched with metals, and galactic outflow is thought to be 
the r esponsible mechanism to spread me tals into the IGM 
fe.g.. lAguirre et al.ll20Qll : ICen et al.ll2005l ). Therefore, an ap- 
propriate model of galactic outflow has to be implemented 
in all cosmological simulations of galaxy formation in order 
to reproduce the observed statistics of galaxies and IGM. 
There are two major astrophysical sources that drive 
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the galactic outflow: supernovae (SNe) and supermassive 
black holes (SMBHs). The SN feedback probably plays an 
important role throughout the entire history of the universe 
as long as the star formation is ongoing, and it is considered 
to be responsible for enriching the IGM. The SMBH feed- 
back is considered to play an important role in determining 
the e volution of massive galaxies at low-z fe.g., [Keres et al. 
l2009l ). In this paper, we focus on the outflows driven by 
SNe, and examine the effects of different models of galactic 
outflows on galaxy growth and IGM enrichment. 

Recent observations show that the galactic outflow 
consists of multiphase gas. Soft X-ray observations reveal 
that some galaxies are sur rounded by haloes of hot gas 
([Strickland Stevensllioooh . A diffuse hard X-ray emission 
that would be associated with hot plasma of temperature 
T = 10^ — 10^ K was detected hy Chan d ra X-ray Observator 
in M82 ([Strickland Heckmanll2Q09l ). iMartiJ (l2QQ5l . l200i ' 
studied the Nai absorption lines in ultra- luminous infrared 
galaxies, and argued that a significant fraction of galactic 
outfiow contains relatively cool (or warm neutral/ionised) 
gas. Although there are many ongoing observational stud- 
ies, the detailed composition of multiphase outflow is still 
poorly constrained. 
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The absorption line studies of outflowing gas from star- 
burst galaxies suggest that the outflow velocit y is a function 
of galaxy mass or star formation rate (SFR) ([Martini [20051 : 
iRupke et al.ll2005h . Using a large sample of low-resolution 
spectroscopy for z ^ 1 galaxies from the DEEP 2 survey, 
IWeiner et all (l2009l ) argued that the outflow velocities of 
z 1 galaxies scale as Vwind oc SFR°"^. iMartinI (|2005h also 
reported that t^wind oc SFR°'^^ These observations attempt 
to make connections with theoretical mo dels of wind driv- 
ing mechanisms fe.g., iMurrav et aI]|2005l V but a clear con- 
nection between the theory and observations is yet to be 
established. 

Our current understanding of the wind driving mech- 
anisms and the propagation processes of outflow from first 
principles is still crude. Although there have been some re- 
cent progress using high-resolution hydrodynamic simula- 
tions in generating galactic outfl ows and res olving the evo- 
lution of shocked, swept-up s hell (iMac Low Ferrara 1993; 
iMori et al]|2002l : iFuiita et al.ll2009l V manv uncertainties re- 
main regarding the energy source of the wind, instabilities 
in the outflowing medium, and the interaction between the 
wind and ambient gas. 

Moreover, in cosmological hydrodynamic simulations, 
there is a so-called "overcooling" problem: if the SN feedback 
energy is deposited as a thermal energy and the simulation 
doesn't have an adequate resolution to resolve the expand- 
ing hot bubble of gas, the thermal energy is quickly radiated 
away without having any effect on the surround ing medium 
re.g.. lKatz eraIlll996l : [Ceverino Klvpinll2009l ) . In the real 
Universe, overlapping SN explosions form a low- density su- 
perbubble of hot gas, which lasts longer and significantly 
affects the evolution of surrounding medium by suppressing 
subsequent cooling of gas. These hot superbubbles cannot 
be resolved even in the highest resolution cosmological sim- 
ulation of today's, if one wants to follow the dynamical evo- 
lution of numerous galaxies simultaneously in a comoving 
volume of >10Mpc. Therefore, the present-day cosmologi- 
cal hydro simulations could still suffer from the overcooling 
problem if no further treatment is given. 

To avoid this problem, many cosmological simula- 
tions adopt a phenomenological model of galactic outflow. 
One widely adopted treatment employs a subgrid multi- 
phase ISM model and the ejection of SPH particles as 
winds (ISpringel HernauistI (20031 : iTornatore et al.1 l2004l : 
lOppenheimer Davell200i )7 In this model, the SN energy is 
deposited in the hot ISM, and the galactic outflow is gener- 
ated by temporally converting some gas particles into wind 
particles based on the mass-loading factor: 



(1) 



where Mw and are the wind mass transfer rate and 
galaxy SFR, respectively. The wind particles get the ini- 
tial kick according to the amount of SN energy injection, 
and all wind particles initially have the same velocity and 
the same 77, i.e., the "constant velocity wind" model. During 
the wind phase, the particles are temporally decoupled from 
the hydro force and are able to escape from the galaxy. This 
wind particle method can enrich the IGM, and regulate the 
galaxy growth and star formation. 

We point out that there are several other attempts 
to model galaxy outflows in cosmological simulations. 
IStinson et al.l (|2006l ) proposed to suppress gas cooling during 



the expansion of SN bubble phase, whose size is computed by 
a blast- wave model. Thi s approach allows them to avoid the 
"overcooling" problem. IScannapieco et aP (|2006l l proposed 
to explicitly divide the SPH particles into cold and hot phase 
particles, and performed the SPH smoothing separately for 
each phases to avoid the artificial cooling caused by the 
mixing of different phases. Their model also suppresses the 
cooling of SN energy injected into the cold phase, and con- 
verts the cold phase to t he hot phase after a cert a in am ount 
of SN energy injection. iDalla Vecchia Schavd (|2008h im- 
plemented the kinetic feedback by injecting the SN kinetic 
energy into the neighbouring gas particles. Each of these 
attempts have remedied the overcooling problem to some 
extent and improved the galactic outflow model in cosmo- 
logical hydrodynamic simulations. 

Although the w i nd pa rticle method developed by 
ISpringel HernquistI (|2003l ) and other alternative ap- 
proaches have been successful to some level, there are still re- 
maining issues. First, there are two theoretical ideas for the 
wind driving mechanism — the "energy- driven" wind and 
the "momentum- driven" wind model. In the energy-driven 
model, the thermal energy injected by SNe into the ambient 
medium drives the galactic wind through the expansion of 
hot bubbles. In the momentum-driven model, the momen- 
tum applied by the radiation pressure from massive stars 
pushes the dus t particles, which e ntrain the gas to drive 
galactic winds (Murray et al. 2005). An important point is 
that the two models have different values of 77 parameter: 



rj oc 



(energy — driven) 
(momentum — driven). 



(2) 



where cTgai is the velocity dispersion of the galaxy 
(iMurrav et al.ll2005f ). However, we note that these relation- 
ships are derived from simple scaling arguments and may not 
be so robust. The values of 77 ar e only loosely constrained by 
observations (e.g.. lMartinlll999l l. The relative importance of 
the two theoretical ideas has not been considered carefully 
in cosmological hydrodynamic simulations. 

Second, there are many problems in the previously 
adopted wind models. For example, the constant velocity 
wind model seems to predict too little C iv absorption, too 
small Clll/Civ ratios ([Aguirre et al.ll2005h , and overheat the 
IGM (jOppenheimer Sz Davell2006h . The latter work also ar- 
gued that the momentum-driven wind model produces more 
mass-loaded but low velocity winds from early galaxies, so 
as to eject a substantial metal mass into the IGM without 
over- heating it. 

Third, the multiphase nature of the wind medium has 
not been recognised very much for the galactic wind models 
in cosmological simulations. In the earlier work, the wind 
particles were selected only from star-forming particles or 
their neighbours. Owing to the hot gas heated by the SNe in 
the subgrid ISM model, the gas temperature in and around 
the star- forming region is very hot. In addition, the high 
constant wind velocity may have heated the IGM too much, 
and we need to revise the wind model to incorporate a cold 
component into the wind medium. 

The best way to resolve these issues is to perform hydro- 
dynamic simulations with a very high-resolution, and simu- 
late the galactic outflow in a cosmological environment from 
first principles. However, the computational expense of such 
a cosmological simulation is still very high and not possible 
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at this moment. An alternative approach is to develop more 
realistic phenomenological galactic outflow models, verify 
the model by comparing with observations, and we can gain 
deeper insight into the phenomenon of galactic outflow. In 
this paper, we adopt the latter approach to study the role of 
galactic outflow in galaxy formation, and develop a new phe- 
nomenological wind model. The new outflow model relies on 
the observational constraints: the wind velocity (vwind) as a 
function of galaxy SFR, and 77 based on galaxy mass. Using 
this phenomenological approach, the new model captures 
the multicomponent origin of wind material and variable 
wind velocity. 

This paper is organized as follows. In Section [2l we de- 
scribe our simulation with a focus on the new wind model. 
In Section [3l we present some of the global statistics such 
as the overall distribution of gas temperature and metals, 
cosmic SFR, and galaxy mass functions. In Section [J we 
focus on the IGM properties such as the mass fractions in 
various phases, overdensity-metallicity relationship, and the 
evolution of cosmic CiV mass density. In Section [5l we will 
summarise our findings and discuss the effects of our new 
wind model. 



2 NUMERICAL TECHNIQUE 
2.1 Simulation Setup 

We will use the modified version of the Tree-particle-mesh 
(TreePM) smoothed particle hydro dynamics (SP H) code 
GADGET-3 (originally described in lSpringelllioOSI ). In this 
code, the SPH calculation is performed based on the en- 
tropy conservative formulation (Springel & Hernquist 2002). 
Our conventi onal code includes radiative cooling by H, He, 
and metals (|Choi k, Nagamin3 l2009bl). heating by a uni - 
form UV b ackground of a modified iHaardt Madaul (|l996l ) 
spectrum (|Katz et al.1 Il996l : IPave et al.1 Il999l ) , star forma- 
tion, supernova feedback, phenomenological model for galac- 
tic winds, and a sub-resolu tion model of multiphase ISM 
(ISpringel k Hernauistll2003l ). 

In the multiphase ISM model, high-density ISM is pic- 
tured to be a two-phase fluid consisting of cold clouds in 
pressure equilibrium with a hot ambient phase. Cold clouds 
grow by radiative cooling out of the hot medium, and this 
material forms the reservoir of baryons available for star for- 
mation. 

For the star formati on model, we adopt the "Pres- 
sure model" described in IChoi Nagamind (|2009al ). This 
star formation model estimates the SFR based on the 
gas pressure rather than gas density, and implicitly con- 
siders the effect of H2 for n iation on star formatio n 
(jSchave Dalla Vecchial I2OO8I: IChoi Nagami ne 2009a). 
Interested readers can refer to lChoi Nagamina (;2009b .a) 
for detailed descriptions. 



2.2 Multicomponent Variable Velocity Wind 
Model 

We will call our phenomenological model of galactic wind 
as the "Multicomponent Variable Velocity" (MVV) wind, 
because the wind material is taken from both low tem- 
perature, high density regions and high temperature, low 



density regions. The fundamental parameters of our MVV 
model model are 77 and t'wind, which are constrained by the 
observations. As mentioned in Section [T] there are numer- 
ical limitations as to properly reproduce the formation of 
galactic outflow. In addition, newly proposed radiation pres- 
sure model is based on a significantly different wind launch- 
ing mechanism compared to the well-developed supernova- 
driven wind. Therefore, we use observational constraints to 
compute Tj and Vwind- 

Most observations provide the values of 77 and Vwind 
as a function of host galaxy mass/SFR. To compute the 
galaxy mass and SFR for galactic winds, we have imple- 
mented a group-finder algorithm into GADGET-3 to com- 
pute the galaxy masses and SFR on-the-fly while the simula- 
tion is running. The group-finder, which is a simplified vari- 
ant of the SUBFIND algorithm developed bv lSpringel et al.l 
(|200lh . identifies the isolated groups of star and gas parti- 
cles (i.e., galaxies) based on the baryonic density field. The 
detailed procedure of this galaxy grouping is described in 
iNagamine et al.l (|2004l V The outer baryonic density thresh- 
old for a galaxy is O.Olpth, where pth is the threshold density 
for star formation. The group properties such as the stellar 
mass, SFR, and maximum size of the galaxy are distributed 
to and stored by the individual member particles. During the 
course of simulation, the group-finder is invoked whenever 
the total SFR or total stellar mass increases by 10% from the 
previous group finding. In addition, the group finding is also 
performed when the snapshot file is generated. In particular, 
the frequency of group finding (Az) is approximately 0.11 
for the '1.5ME' run in the redshift range 3 < z < 5. The 
wind parameters (77 and t'wind) are computed based on the 
group properties stored in the member particle data, as we 
describe below. 

In our MVV wind model, the chance that a gas particle 
becomes a wind particle is based on the values of 77, galaxy 
mass, and galaxy SFR. The model assumes that all gas par- 
ticles in a given galaxy have the same chance to become part 
of the wind, which has an important consequence. Since any 
gas particle inside a galaxy, either hot or cold, could become 
a wind particle, the galactic outflow naturally reproduces a 
multiphase nature of outflow medium. The multiphase na- 
ture of galact ic outflow is suggested by many o bservat ions 
(lMartinll2005l . I2OO6I : IStrickland Heckmanll 20091 ) , but it has 
not been taken into account explicitly in cosmological hydro- 
dynamic simulations. In our model, the wind particles are 
randomly chosen from constituent gas particles of a galaxy, 
because we currently have very limited knowledge on the 
origin of the multiphase outflow medium. The outflow in 
our model includes gas from high density, cold gas in the 
galaxy, resembling the cold component entrained by the hot 
expanding gas from SNe fe.g.. tHeckman 2002i V Note that the 
'multiphase' nature of the MVV wind model does not mean 
that a single wind particle is made of two or three phase 
medium, but it simply means that the new wind model gen- 
erates a wind made of both hot and cold gas. 

In the MVV wind model, Vwind depends on the galaxy 
mass and galaxy SFR. We flrst assume that the wind veloc- 
ity is proportional to the escape velocity of the galaxy: 



'^wind 



(3) 



where ^ is a scaling factor of wind velocity. We then employ 
an empirical relationship between Ve^c and galaxy SFR as 
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Figure 1. Input wind parameters ('Uwind cind 77) for the MVV wind model. Left: Wind velocity as a function of galaxy SFR assuming 
C = 1. We compare the -Uwind at z =3 and 5, as well as the constant wind velocity (484 km s"-*^) in the CW run. The yellow shade is the 
observed wind velocity range from lWeiner et al.l (l2009h . Right: Wind mass-loading factor as a function of galaxy velocity dispersion. We 
compare 77 for the energy driven wind and momentum driven wind, as well as the constant 77 in the CW run. The points are the observed 
data from I Martini (| 19991 ). 



follows. From the halo mass (Mhaio) 
circular velocity as 



we can compute the 



= ( 



\ R200 



1/2 



124 



Mha 



10ii/i-iM( 



1/3 



2/3 



47r _ 

-^"^200 



1/3" 



1/2 



1/2 



km s 



(4) 



where R200 and p2oo is the radius and mean density of the 
halo with an overdensity of 200. 

The SFR can be estimated from the UV l uminosity: 
SFR[MQyT-^] = 1.4 X 10"^^L^ [ergs"^ Hz"^] (iKennicuttI 
1 19981 ). Earlier, using the res ults from cosmological simula- 
tions, iNagamine et al.l (|2007l l obtained the relationship be- 
tween the absolute AB magnitude of a galaxy and halo mass: 
Mab = —2. 5 (log Mhaio — 12) — 21.5, where Mhaio is in units of 
h~^MQ. Combining the above two relationships, we obtain 



\og{SFR) = log(Mh 



10.62. 



(5) 



This relationship provides a crude, but a reasonable conver- 
sion between Mhaio and SFR. 

Combining Equations and (O, and using Vesc = 
V^Vc, we obtain the relationship between Vesc and SFR as 



130 (SFR) 



1/3 



1/2 



km s 



(6) 



The above relatio nship Vesc oc (SFRY^^ is consist ent with 
the observed one (|Martinl 120051 : 1 Weiner et al.ll2009l ). 

As discussed in Section [U the two wind-driving mech- 
anisms suggest two different scalings for 77. In the energy- 
driven wind case. 



77 = (cro/cTgal) , 

where ctq = 300 km s~^, and (jgai 



(7) 

c/2 is the velocity 



dispersion of the galaxy. In the momentum- driven wind case, 



r] = cro/cTg 



(8) 



In Figure [T] we show the scaling of 77 and t'wind as a function 
of galaxy SFR and velocity dispersion a. In the left panel, 
one can see that, in the MVV model the wind velocity is an 
increasing function of SFR and redshift, in contrast to the 
constant wind velocity model. Unfortunately, it has not been 
understood which wind-driving mechanism is the dominant 
process, and the current observational data is not yet robust 
enough to clearly discern the best wind model. Therefore, 
we develop a new wind model by employing both mecha- 
nisms in this paper. In the low-density regions (pgas < Pth), 
the gas is sufficiently far away from star-forming regions, 
and the hot bubbles of gas are easy to expand and escape 
the region. Therefore in our fiducial wind model, we em- 
ploy the energy-driven mechanism for the gas particles with 
Pgas < pth- On the other hand, in the high-density star- 
forming regions (yOgas > Pth), gas is close to the massive 
stars and SNe, and more likely to be subject to a strong 
radiation from these sources. Therefore in our fiducial wind 
model, we employ the momentum-driven mechanism for the 
gas particles with pgas > Pth- 

Since the new model generates the wind particle from 
both star-forming and non-starforming regions, the 77 pa- 
rameter can be written as an addition of two separate parts: 



r] = VSF + r]n 



-SF, 



(9) 



where r/sF — M-w,sf / for the wind from star- forming 
gas, and r]non-SF = Mw,non-SF / for the wind from non- 
star-forming gas. Therefore, 



-SF - (VSF + r}non-SF)M^ (10) 



for the total wind. Since the wind particles are selected ran- 
domly within a given galaxy, the mass fraction of the wind 
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High density Low density 



Model 




c 




c 


1.5ME 


Momentum 


1.5 


Energy 


1 


ME 


Momentum 


1 


Energy 


1 


EE 


Energy 


1 


Energy 


1 


1.5MM 


Momentum 


1.5 


Momentum 


1.5 


NW 


No 


galactic wind 





CW Constant velocity wind 



Table 1. The six different wind models employed in this paper. 
The parameter 77 is the mass-loading factor, and (" is the scaling 
parameter for the wind velocity. As for the model name, for ex- 
ample, "1.5ME" consists of the value of "M(omentum)" , and 
"E(nergy)", for high and low density gas, respectively. 



material that comes from star-forming and non-starforming 
gas would depend on the masses of the two components. 

We note that our new wind model is an alternative 
phenomenological model for galactic outflow in cosmolog- 
ical SPH simulations. The values of wind parameters are 
determined by observational constraints rather than from 
physical first principles. In addition, the wind particles are 
decoupled from the hydro force during the wind phase. 
This assumption has been cri ticised by a few authors (e.g., 
iDalla Vecchia Schavdlioosl ) owing to its unphysical na- 
ture. However, it is expected that the hydrodynamical inter- 
action between wind particles and other gas particles in the 
galaxy could give rise to unwanted effects for low-resolution 
simulations. If galactic outflow cannot be launched due to 
hydrodynamic interaction, there is no point in implement- 
ing this phenomenological model in our simulation. In the 
future, we need to continuously improve our model so that 
the wind generation and propagation will be based on phys- 
ical principles. 

In the default model, we adopt C=l-5 for the 
momentum-driven wind (for the high-density gas source) 
and C=l for the energy-driven wind (for the low-density gas 
source). We call this default model the "1.5ME" run. The 
value of C for the momentum-driven wind can be greater 
than unity, because the radiation pressure applied by the 
massive stars can continuously push the outflowing gas to a 
velocity greater than the escape velocity. Besides our fidu- 
cial wind mixture, we also test several other mixtures of 77 
and C, as shown in Table [T] For comparison purpose, we also 
consider the model s with no wind (NW) and a constant ve- 
locity wind (CW) (ISpringel k Hernauistlliooi ). In the CW 
run, only star-forming SPH particles can become wind par- 
ticles, and fixed values of 77 = 2 and t'wind = 484kms~^ are 
adopted. 

In this paper, we adopt the latest WMAP5 cosmology: 
= 0.26, r^A = 0.74, = 0.044 h = 0.72, = 0.96, 
and (78 = 0.80 (|Komatsu et al ] l2009l ). The simulations are 
initialised with 216^ particles each for gas and dark matter 
in a comoving box of (10 Mpc)^. We stop this series at 
z — 2.75, as it misses the long wavelength perturbations at 
lower redshifts. In these simulations, dark matter particle 
mass is 5.96 x lO^/i^^M©, and the initial gas particle mass 
is 1.21 X 10^h~^ Mq. Using the identical initial condition for 
all the runs, we examine the effects of different wind models 
listed in Table [T] 



3 GLOBAL PROPERTIES 

3.1 Distribution of Gas, Thermal Energy, and 
Metals 

Figures [2] - 2] show the snapshot view of the simulations 
with six different wind models at z=?>\ the CW, NW, ME, 
1.5ME, EE, and 1.5MM run from top-left to bottom right. 
Figure [2] shows that the gas column density distribution is 
very similar in all the runs, which implies that the amount 
of outflow material is negligible compared to the total gas 
content in the entire simulation box. 

In contrast, the effect of outflow on the IGM temper- 
ature and the global metal distribution is very significant. 
Figure [3] shows that there is a clear distinction between the 
CW run and the others in the distribution of gas thermal en- 
ergy. In the CW run, the IGM is significantly heated and the 
internal energy is much higher on larger scales than the other 
runs. The similar energy distribution between the NW run 
and the MVV runs suggests that the MVV wind models can 
generate galactic outflows without overheating the IGM, not 
as much as the CW run. This is an encouraging feature for 
the MVV wind model, because lOppenheimer Pavel (l2006h 
showed that the IGM temperature in the CW run is too high 
compared to the observed IGM te mperature measured from 
the line- width of the Ly-a forest (|Schave et al.ll2000l V We 
will discuss this issue quantitatively in Section W72\ 

Figure |4] shows the projected metallicity distribution 
for all the runs at ^ = 3. The CW run spreads the metals 
to the greatest distances, whereas the NW run spreads al- 
most no metals over the intergalactic space. This suggests 
that the mechanisms other than the galactic outflow, such 
as ram pressure stripping, are inefficient in enriching the 
IGM. The MVV runs spread a significant amount of metals 
over the IGM, but not as much as the CW run. These snap- 
shots demonstrate that the MVV wind can enrich the IGM 
without overheating it. 

Among the four MVV wind models, the 1.5ME run 
shows the widest metal distribution. The comparison be- 
tween the 1.5ME and the ME run tells us that the high t'wind 
(C = 1.5) causes a wider metal distribution. The compar- 
isons between the EE, ME, and 1.5MM run demonstrate the 
effect of mass- loading parameter 77: the momentum-driven 
wind is overall more efficient in the IGM metal enrichment 
than the energy-driven wind. The energy-driven wind tends 
to eject more winds from lower mass galaxies, which have 
lower metallicity, than the momentum-driven wind, owing to 
the cTgai dependence of 77. The lower metallicity of the winds 
from lower mass galaxies gives rise to the less metal spread 
in the energy-driven wind than the momentum-driven wind. 
In the MVV wind models, the choices of 77 and t'wind alter 
the degree of IGM metal enrichment. 



3.2 Cosmic Star Formation History 

Figure [5] shows the cosmic star formation history of all the 
runs with different wind models. The MVV wind runs show 
a significant suppression of SFR compared to the NW run, 
but a similar to or a higher SFR than the CW run. Ob- 
servations suggest that the range of co smic SFR at z ^ 3 
is be tween 0.07 - 0.3 M© yr"^ Mpc"^ (|Hopkins Beacom 
l2006h . although there is a huge uncertainty (|Kistler et al.l 



© 0000 RAS, MNRAS 000, 000-000 



6 Choi & Nagamine 




Figure 2. Gas column density at z=3 projected along with 2;-axis for the CW, NW, ME, 1.5ME, EE, and 1.5MM run from the top-left 
to bottom-right. All panels use the same colour scale as shown in the side bar. The size of each panel is comoving 10 /i~-'^Mpc on the 
side. The ga s distribution is almost identical in all the runs. The visualisation was done by the SPLASH code with 3-D rendering mode 
(|Pricdl2007l V 



I2OO9I : IChoi Nagamin3l2QQ9al ). The comparison with the 
observed cosmic SFR shows that the NW run generates too 
many stars, while the results of the MVV wind runs are in 
the observed range. 

The differences in cosmic SFR among the MVV runs 
yield insights on the effects of wind parameters (77 and C)- 
Among the MVV runs, the EE run shows the lowest SFR 
and a good agreement with that of the CW run. As discussed 
above, the energy-driven wind tends to produce galactic out- 
flows more in lower mass galaxies. Since low- mass galax- 
ies form earlier in the hierarchical galaxy formation model, 
the energy-driven wind produces more outflows in the ear- 
lier epoch of galaxy formation. Consequently, the energy- 
driven wind is more efficient in suppressing SFR than the 
momentum-driven wind. The 1.5ME run suppresses the SFR 
more significantly than the ME run, because of the higher 
value of Vwind- 

Interestingly, the cosmic SFR in the MVV runs does 
not show a turn-over at 2: > 3. It implies that the peak 
of the cosmic SFR is located at 2; < 3 for the MVV runs, 
while the CW run has a clear peak at z > 3.5. As shown 
in Figure O the CW model overheats the IGM, which pro- 
hibits the gas accretion onto galaxies and suppresses star 
formation. On the other hand, the MVV wind model does 
not heat the IGM as much as the CW model, and allows 
continuous IGM accretion that increases the cosmic SFR. 
Therefore, the MVV wind model can shift the peak of the 
cosmic SFR to a lower redshift without changing the star 



formation model fsee IChoi Nagamind [ioOQal . for the ef- 
fect of star formation model on the cosmic SFR). In order 
to definitely locate the peak of the cosmic SFR, we need 
to run simulations with a larger volume which reach lower 
z. Our current study suggests that the MVV wind model 
changes the shape and the amplitude of cosmic SFR. 



3.3 Galaxy Stellar and Baryonic Mass Functions 

Figure [6] shows the stellar and baryonic (i.e., star+gas) mass 
functions of the six different wind models at 2; = 5, 4, and 
3. In general, the stellar mass functions of MVV runs are 
located in-between the CW and NW run. The trend is con- 
sistent with the cosmic SFR; the CW run generates the least 
amount of stars, while the NW run generates too many stars. 
The consistency between the cosmic SFR and the galaxy 
mass function also holds among the MVV runs: the number 
of galaxies in the EE run is the lowest, and that of the ME 
run is the highest. As discussed above, the choices of 77 and 
Vv^ind parameters determine the amount of star formation, 
as well as the stellar masses of galaxies. 

Interestingly, the low-mass-end slope of both stellar and 
baryonic mass functions in the MVV runs are flatter than 
the other runs. For example, the power-law slope is a ~ —1.6 
(assuming the Schechter function) in the mass range of 
Mstar - 10^ - lO^M© for the MVV run, while a - -2 
in the CW and NW run at 2; = 3. The number of low-mass 
galaxies is reduced in the MVV runs compared to the CW 
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Figure 3. Same as in Figure [3 but for the projected internal energy (/ udz): CW, NW, ME, 1.5ME, EE, and 1.5MM run from the 
top-left to bottom-right. The IGM temperature is significantly higher on much larger scales in the CW run (top-left) than in the other 
runs. The MVV runs hardly heat the IGM on large scales. 



run, because the values of 77 are higher for lower mass galax- 
ies in the MVV runs, while the CW run adopts a constant 
value of 77. Therefore, the low-mass galaxies in the MVV 
runs lose comparatively more mass than those in the CW 
run (see also Section l3.4p . Furthermore, less IGM heating 
in the MVV runs allows continuous gas accretion onto more 
massive galaxies, and the galaxies with Mstar > 10^ M© in 
the MVV run become more massive than those in the CW 



3.4 Gas Fraction of Galaxies 

Figure [71 shows the gas fraction of galaxies at 2; = 5 and 
3. In general, the gas fraction in the MVV run is roughly 
constant as a function of galaxy mass for galaxies with = 
10^ — 1O^M0. At z = 5, the gas fractions of low-mass galaxies 
in the MVV runs are lower than those of the CW and NW 
run, and higher for massive galaxies. Several mechanisms 
are at work to give rise to this result, but the dependency 
of r] on galaxy mass plays a key role, because it is directly 
related to the gas removal rate. Compared to the CW run, 
the MVV runs have higher 77 and more gas is ejected from 
lower mass galaxies. Compared to the NW run, the MVV 
runs have lower SFRs, therefore they consume less gas in 
massive galaxies, resulting in higher /gas in massive galaxies. 
Among the MVV runs, the EE run shows the lowest gas 
fraction for the lower mass galaxies and highest gas fraction 
for the higher mass galaxies, because 77 is the greatest in 
the EE run for low- mass galaxies and smallest for massive 



galaxies among different MVV runs. Note that not only 77 
but also fwind characterise galactic outflows. However, t'wind 
mostly affects the IGM properties. Therefore, we argue that 
the difference in /gas mainly comes from the 77, which is more 
closely related to the gas removal rate. 

At z — 3, the general trend is similar to z — 5. How- 
ever, at lower redshifts, the balance between outflow and 
IGM accretion become more important for the gas fraction 
distribution. At z = 3, the CW run has the lowest /gas, be- 
cause its high IGM temperature reduces the IGM accretion 
onto galaxies. One interesting feature in the z = 3 result is 
that the EE run shows the highest /gas- This is because the 
SFR is the lowest in the EE run as we saw in Figure [5l and 
more gas is preserved in galaxies. The IGM temperature is 
also lower than the CW run, enhancing the IGM accretion 
onto galaxies. We will also see the enhanced IGM accretion 
in the EE run compared to the CW run in Section 14.11 



3.5 Wind velocity 

Figure [8] shows the wind velocity as a function of galaxy 
SFR. We directly measure the velocities of the wind particles 
for each galaxy, and then average them over for each galaxy. 
We assume that the wind particles were ejected from the 
nearest galaxy. This assumption may lead to some wrong 
associations between the wind particles and the launching 
galaxy, but the fraction of such mis-identification must be 
small and should not affect the average velocity. The figure 
shows that the wind velocity of MVV runs is a function of 
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Figure 4. Same as in FigureE] but for the projected metallicity (/ Zdz): CW, NW, ME, 1.5ME, EE, and 1.5MM run from the top-left 
to bottom-right. The metals are spread over much broader regions in the CW run than in the other runs. The MVV runs also spread a 
significant amount of metals over the IGM, while the NW run (top middle panel) hardly spreads metals. 
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Figure 5. Cosmic star formation history of the simulations with 
six different wind models. The MVV wind runs have significantly 
lower SFRs compared to the NW run. The CW and EE run have 
the lowest cosmic SFR. For comparison, the yellow shaded re- 
gion indicates th e locu s of the observed data points compiled by 
iNagamine et aP (l2006l V 



SFR, in good agreement with the input velocity (Eqns. [3l[6]; 
indicated by the long-dashed line in the top left panel), while 
the wind velocity of the CW run is independent of galaxy 
SFR. Both the MVV and CW runs show a non-negligible 
scatter, which results from the velocity dispersion of gas 
particles in a galaxy prior to becoming a wind. 

W e compare our results with the observed wind velocity 
bv lWe iner et al. (2009). Owing to the flux limit, the observa- 
tion only probes the wind velocity of massive galaxies with 
SFR greater than 10 Mq yr~^. The observations of wind ve- 
locities also contain significant scatter , and it is difficult 
to distinguish a favored model ([Martini [2005"; 'Wei ner et al.l 
[2009). However, many observations have suggested that the 
wind velocity is a function of galaxy size with a similar re- 
lationship to our model, and the results of the MVV run is 
consistent with the limited observations of wind velocities. 

The constant wind velocity in the CW run causes the 
high IGM temperature (see Figure [3] and Section [4. 2p . How- 
ever, this difference in the IGM temperature does not result 
from different outflow energy. The injected kinetic energies 
depend on both the outflow velocity and the amount of the 
outflow mass, which depends on 77. We find that the injected 
kinetic energies by winds in the CW run and the MVV 
runs for a given galaxy SFR are indifferent, because low- 
mass galaxies in the CW run generate high- velocity outflow, 
while the low-mass galaxies in the MVV runs have higher 
mass loading factors. What makes a big difference in the 
IGM temperature is the high wind velocity from low-mass 
galaxies. Since the potential wells of these galaxies are shal- 
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Figure 6. Evolution of galaxy mass functions for the simulations with six different wind models. The left panels are the stellar mass 
functions, and the right panels are the baryonic (star + gas) mass functions. The mass functions in the MVV runs are flatter than that 
in the CW and NW runs. 



low, their winds propagate to fill a significant volume, and 
their effect on the environment is significant. In contrast, the 
winds from low-mass galaxies in the MVV runs do not have 
such a high velocity, so their effect on the environment is 
limited. The high IGM temperature in the CW run conflicts 
with the observed IGM temperature. 



4 EVOLUTION OF IGM 
4.1 Phase Evolution 

Figure [9] shows the evolution of mass fractions of different 
baryon phases. The baryons in the Universe can be broadly 
categorised into four different phases according to their 
overdensity and temp erature : ^hot ', Varm-hot', ^d i ffuse' , 
and ^condensed' (Dave et al.1 Il999l : ICen k Ostrike3 Il999l : 
iDave et aP l2QQlh . The 'hot' phase is the non- star forming 
gas with T > 10^ K. The 'warm-hot' phase is the non- 
starforming gas with 10^ < T < lO'^ K. The 'diffuse' phase 
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is non-starforming gas with T < 10^ K and S < 1000, where 
6 = p/ p — 1 and p is the mean baryonic density. The 'con- 
densed' phase includes non-starforming gas with T < 10^ K 
and 6 > 1000 , as well as the star-forming gas and stars. The 
'hot' and 'warm-hot' phases are mostly the shock- heated 
gas in clusters and groups of galaxies. The 'diffuse' phase is 
mostly the photoionised IGM with lower temperature, which 
can be observed as the Ly-a forest. The 'condensed' phase 
is the high density, colder gas in galaxies. 

The panels (a) and (h) show that the MVV runs always 
contain more diffuse gas and less 'warm- hot + hot' gas than 
the CW run. This difference results from the velocity and the 
temperature of the outflowing gas. After the wind particle is 
ejected, it is mixed with the ambient gas and the tempera- 
ture can be decreased via adiabatic cooling. Meanwhile, the 
ambient gas temperature increases due to the thermal en- 
ergy from the outflowing gas and shock heating. As discussed 
in Section [331 the galactic outflow from low-mass galaxies 
in the CW run influences the ambient gas significantly. In 
general, the heating from these outflows increase 'warm-hot 
+ hot' gas fraction in the CW run. The comparison of the 
CW and MVV run in panels (a) and (h) shows that most of 
the outflow medium goes into the 'warm- hot + hot' phase 
in the CW run, whereas in the MVV runs, it goes into the 
'diffuse' phase. In conclusion, the MVV wind model signifi- 
cantly alters the mass fraction of different baryon phases. 

The evolution of the 'condensed' and 'IGM' phase {pan- 
els (c) and (d)) is related to the cosmic SFR, because the 
amount of 'condensed' phase and SFR is closely correlated. 
The CW run shows the highest fraction of the IGM and 
the lowest fraction of the condensed gas, while the NW run 
shows the opposite. Among the MVV runs, the EE run has 
the highest fraction of the IGM, owing to the stronger out- 



flows from low- mass galaxies. We also see that the IGM mass 
fraction rapidly declines in the EE run at z < 4 compared 
to the CW run. This is because the low SFR produces weak 
galactic outflow and the low IGM temperature increases the 
IGM accretion onto galaxies in the EE run. This increase of 
the IGM accretion is also responsible for the high /gas for 
the galaxies in the EE run at z = 3 (see Section [X4l 

4.2 Temperature Evolution 

Earlier, we showed in Figure [3l that the NW and MVV run 
hardly heat the IGM, while the IGM in the CW run is signif- 
icantly hotter. In order to quantitatively discuss this issue, 
we plot the evolution of volume- weighted temperature of all 
gas in the simulation box in Figure [TOl Note that the overall 
IGM property is better represented by the volume- weighted 
value, because the IGM covers a large volume of the Universe 
including low-density voids and high-density filaments. Fig- 
ure (TOl confirms that the IGM temperature in the CW run is 
higher than the other runs at all redshifts. It confirms that 
the average IGM temperature in the MVV run is very sim- 
ilar to that of the NW run. Since the IGM heating by the 
outflows in the NW run should be minimal, this figure also 
shows that the MVV wind model indeed does not heat the 

IGM^ 

ISchave et al.l (|2000l ) estimated the IGM temperature at 
the mean density using the Ly-a absorption line widths, as 
shown by the yellow shade in Figure [TOl At 3.5 < z < 4.5, 
the results of the MVV runs are at the lower edge of the 
observed estimates. At z < 3.5, the observed temperature 
range of ^ 1.5 — 2.5 x 10^ K is higher than the MVV run 
results 1 X 10^ K), and similar to or lower than the CW 
run result (- 3 x 10^ K). ISchave et all tOOd ) argued that 
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Figure 8. The averaged outflow speed as a function of galaxy SFR. For each galaxy, we measure the wind velocity of all wind particles 
and compute the average wind velocity, as represented by the red dots. The blue solid lines are the average speed for galaxies within a 
given mass range. The black dashed line is the input wind velocity wit h C = 1 at z =3, whi ch is identical to the z = line on left panel 
of Figure [T] The yellow shade is the observed wind velocity range from lWeiner et a" 



the IGM temperature shows a peak and the gas becomes 
nearly isothermal at 2; 3, which can be interpreted as 
an evidence for the Hell r eionization (e.g., ISokasian et al.l 
I2OO3I : iMcQuinn et aLllioOol ). The optically thin approxima- 
tion is adopted for the UV background radiation in many 
cosmological simulations, however, this method underesti- 
mates the Hell photoheating rate. The mean distance of 
He Lyman limit systems is smaller than the mean distance 
between the sources of ionising photons, in which case the 
optically thin approximation is not appropriate, and a sig- 
nificant fraction of emitted ionising photon is absorbed and 
enhances the photoheating rate. We need to consider the 
full radiative transfer effect in order to properly estimate 
the thermal history of IGM around He 11 reionization epoch 
([Abel Haeh nelt 1999|). The neglect of this radiative trans- 
fer effects results in the IGM temperatures which are too low 
by a factor of 2 after Hell reionization (jjena et al ] l2QQ5l : 
iFaucher-Giguere et aLllioOOI V Therefore the IGM tempera- 
ture 2; 3 in simulations with an optically thin approxi- 
mation should not be higher than the observations, if the 
assumed amplitude of UV background is reasonable. Our 
results suggest that the CW run clearly overheats the IGM, 
while the MVV runs do not. 

We also note that the IGM temperatures in the 
CW and NW run are l ower than that reported by 
iQppenheimer k Pavel (1 20061 ). despite their outflow prescrip- 
tion is very similar to ours. This may result from the different 
cosmology such as lower as and lower ris that we adopt, or 
diflPere nt metal cooUn^ m ethod from theirs. In particular, re- 
ce ntlvlShen et al.l (l2009l) argu ed that the metal cooling rate 
in IChoi k, Nagamind ( 2009bl l was overestimated compared 



to their more reflned treatment, owing to the missing treat- 
ment of metal ionisation by the UV background effect. We 
will investigate this issue in the future. 

4.3 Metal Enrichment 

4.4 Overdensity — Metallicity Relationship 

Metals are generated in galaxies that are usually located in 
dense regions, and spread over the low-density IGM. There- 
fore it is expected that there will be a positive correlation 
between gas density an d metallicity. Using the C IV pixel 
optical depth statistics, ISchave et al.1 (I2OO3I ) estimated the 
relationship: [C/H] = -3, 47+0.08(^-3) + 0.65[log(5) -0.5], 
where z is the redshift range between 1.8 to 4.1, and 5 = 
p/p— 1 is the baryonic overdensity. This relationship indi- 
cates a positive correlation between the IGM metallicity and 
baryonic overdensity. 

Figure [TT] shows the baryonic overdensity— metallicity 
relation at 2; = 3 for the simulations with six different 
wind models . We a lso overplot the relationship derived by 
^Scha ve et al.l (|2003l l with its 1-cr scatter (shown by the yellow 
shade). The runs with galactic outflows overpredict the gas 
metallicity compared to the observational estimate, while 
the NW run underestimates the metallicity. This suggests 
that the galactic outflow is necessary to reproduce the ob- 
served IGM metallicity, but also that the current outflow 
models are spreading the metals too efficiently. 

In Figure fTTl the CW and NW run are definitely outside 
the 1-cr range of the observed data. The MVV runs gener- 
ally show somewhat higher metallicity than the observed 
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Figure 9. Evolution of the mass fractions of the four different baryon phases (hot, warm-hot, diffuse, and condensed) for the simulations 
with six different wind models. Panel (h) shows the addition of the two phases, 'warm-hot + hot', and panel (c) shows the total mass 
in the IGM, i.e., 'diffuse + warm-hot + hot'. This figure shows that the outflow medium in the MVV wind model mostly goes into the 
'diffuse' phase instead of the 'warm-hot' phase, whereas in the CW run, it is the opposite. 



data, but in better agreement with the data than the CW 
and NW run. It implies that the MVV wind reduces the de- 
gree of metal spreading, more consistently with the observa- 
tion than the CW or NW run. Among different MVV runs, 
the EE run shows better agreement with the observation, 
because in this run, the low-mass, low-metallicity galaxies 
have high mass- loading factors for the energy-driven wind, 
and the outflow medium tends to have a lower metallicity. In 
addition, the cosmic SFR in the EE run is the lowest among 
different MVV runs. The low SFR results in less metal pro- 
duction, which also contributes to the better agreement with 
the observed density-met allicity relationship for the EE run. 



4.5 C IV Statistics 

The most commonly encountered metal lines associated with 
the Ly-g forest are the Civ AA 1548.204 1, 1550.7812 dou- 
blet (ICowie et al.lll995l : lEllison et al.ll2000l ). Therefore, many 
IGM metallicity studies use the Civ measurement s to es - 
timate the total meta ll icity including Sch ave et al.l (|2003h . 
lOppenheimer Pavel (|2006l ) argued that CiV is not a 
straightforward measure of the IGM metallicity owing to 
the evolution of ionisation fraction, which is affected by the 



IGM temperature. They also showed that the observed flat 
evolution of Q(C iv) is caused by the evolution of ionisa- 
tion fraction and the UV background. Although Q.{C) con- 
tinuously increases, the increasing ionisation fraction com- 
pensates, and give rise to a flat evolution of Q(Civ) evo- 
lution. Therefore, comparing the observed IGM metallic- 
ity estimated from the C IV measurements with the total 
metallicity in the simulation may not be a fair comparison, 
and this could be one of the reasons for the discrepancy be- 
tween the observation and simulations shown in Figure [TT1 
To avoid this complication, we directly compare the evolu- 
tion of Q(C iv) in our simulations with observations. 

Owing to the computational cost, our current simula- 
tions do not track the evolution of individual metal ele- 
ment. However we track the evolution of the total metal- 
licity, Z, using the p r imordial metal ratio for Z < 0.1 Zq 
(IWheeler et al.l Il989l : iBessell et all Il99lh, and the solar 
abundance for Z > Zq ([Anders Grevess 3 ll989l ). For 0.1 < 
Z/Zq < 1.0, the abundance pattern is computed by interpo- 
lating bet ween the primordial and t he solar abundance pat- 
terns fsee lChoi Nagamine|[2Q09bl , for detailed discussion). 
We estimated Q(C) from the total metallicity evolution with 
a given abundance pattern. Abundance of C iv is estimated 
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from the abundance of carbon a nd its ionisation frac tion. We 
use the c96 version of Cloudy (|Ferland et alJllQQSl l assum- 
ing an optically thin slab of ^as with the ionisation radiation 
field of iHaardt Madaul (|l996l ). The amplitude of the ioni- 
sation background is reduced by a factor of 1.6 to match th e 
observed Ly-a flux decrement (jOppenheimer Dave|[2QQ6l V 
In Figure [T2l we show the redshift evolution of Q(C) and 
Q(C iv), together with the observed range in yellow shade. In 
our simulations, Q(C) monotonically increases with decreas- 
ing redshift, as the total metallicity increases. The Q(Civ) 
also continuously increases, except that the EE run turns 
over at z < 3.3. The MVV runs show a more mild evolution, 
and Q(C iv) becomes flat at z < 3.5. Observations show that 
the redshift evolution oiQ(C iv) is almost flat between z = 2 
and 4.5, and decreases at 2; > 5 (jRvan- Weber et aP l2009l ) . 
Although the agreement is not perfect, the Q(C iv) evolution 
in the simulations with galactic outflows is more or less con- 
sistent with the observed range (yellow shade) . The NW run 
completely underestimates the Q(C iv) by about an order of 
magnitude, therefore its result is not shown in the lower part 
of Figure 1121 This suggests that a galactic outflow model of 
some form is required to reproduce the observed r2(Civ). 
The flattening of Q(Clv) ai z < 3.5 in the MVV runs is 
encouraging, and these runs may show a better agreement 
with future observations, if we extend the simulations to a 
lower z. 

The discrepancies between our simulations and observa- 
tions in Figure[T2lis significantly less than those in FigurefTTI 
It suggests that the IGM temperature distribution might 
be responsible for the discrepancies in the overdensity- 
metallicity relation as we described earlier. These compar- 
isons suggest that we need to test and improve the UV back- 
ground model as well as the outflow model to resolve the 
discrepancies shown in Figures [TT] and Figure [T2l 

The evolution of Q(C iv) has been used by several previ- 
ous studies to ve rify galactic outflow models in cosmological 
simulations fe.g. JOppenheimer Davi200€ V Compared to 
the Figure 12 of lOppenheimer Pavel (|2QQ6l ). the evolution 
of Q(Civ) in our MVV runs show a similar level of agree- 
ment with observation. Therefore our MVV wind model is 
also a legitimate new approach for modeling galactic outflow 
in cosmological simulations. 



5 SUMMARY AND DISCUSSION 

In this paper, we developed a phenomenological multicom- 
ponent variable velocity (MVV) galactic outflow model for 
cosmological SPH simulations. In this new wind model, the 
wind parameters such as 77 and t'wind, are based on the 
galaxy mass and SFR, which are computed from the on- 
the-fly group-finder. 

In the MVV wind model, all gas particles in a given 
galaxy have the same probability to become winds, and we 
allow lower density, cold gas to be part of the wind, which 
was not allowed in the previous models. This allows our 
MVV wind to reflect the fact that wind material can arise 
from different phases in the ISM, and the MVV wind sig- 
nificantly enriches the IGM metallicity without overheating 
it, unlike the previous CW model. The comparison of our 
simulations with the observed overdensity-metallicity rela- 
tionship shows that the IGM metallicity in our simulation 
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Figure 10. Evolution of the volume-weighted temperature for 
six simulations with different wind models. Except the CW run, 
all simulations show almost identical temperature evolution. The 
yel low shade is t he ra n^e of observed IGM temperature derived 
bv lSchave et al ] (I2OOOI ) from the Ly-o; forest line- widths. 
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Figure 11. Baryonic overdensity-metallicity relation at z=?> for 
six simulations with different wind mo dels. The black solid line is 
the observational estimate derived by ISchave et al.l (I2OO3I ) from 
the Ly-a forest line- widths, and the yellow shade is the 1-cr scat- 
ter. 
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Figure 12. Evolution of (solid lines) and (dashed 

lines) as a function of redshift for all simulations with six dif- 
ferent wind models. The yellow shade ind icates the range of 
observed data of 17(C i v) (fSongaila i200ll : IPettini et all l2003l : 
iRvan- Weber et aljlioool ). Our result of 17(Clv) in the NW run 
is so low that it is located outside of the axis range. 



is slightly higher than the observations. However, our simu- 
lations show a reasonable agreement with the observational 
estimates of Q(Clv), though not perfect. Because current 
observations rely only on a few metal absorption lines such 
as C IV, the effects of UV background radiation and the IGM 
temperature may be responsible for this discrepancy, and we 
need to improve the treatment of UV background in the fu- 
ture. 

Recent WM AP measurement of the electron scatter- 
ing optical de pth (iDunklev et al.l l2QQ9^) and the Ly-a forest 
transmission (|Chiu et al. 12003) suggest that t h e UV back- 
ground radiation model of iHaardt Madau (Il996l) may 
need to be improved. In particular, iFaucher-Giguere et al.l 
(l2009l ) investigated the implications of Ly-a forest opacity 
measurements at 2 < z < 4, and foun d a remarkably f i at ion - 
isation rate over this redshift range. iHambrick et al.1 (I2OO9I ) 
examined the effects of different ionisation rate evolution as 
a function of redshift on galaxy formation and evolution. 
They found that increase in either the intensity or hardness 
of ionising radiation, which recent input UV background ra- 
diation models predict comparing to old models, generally 
reduces early star formation and push it toward lower red- 
shifts. This change also results in enhanced late gas inflows, 
decrease of stellar half-mass radii, increase of central velocity 
dispersion, and reduction of substructures. 

The MVV winds effectively remove the material from 
galaxies into the IGM. In the runs with the MVV wind 
model, the cosmic SFR is significantly reduced compared 
to the NW run, and the results are more consistent with the 
observed cosmic SFR. Because of the variable wind velocity 
and the functional dependence of r\ on cTgai , the galaxy stellar 
mass functions in the MVV runs are flatter than that in the 



CW and NW runs. Owing to the small volume of the sim- 
ulations in this paper, we are only able to demonstrate this 
effect for the high- 2; , low-mass galaxies, and the flux limit of 
current observations do not reach this faintness yet. Never- 
theless, the flat galaxy stellar mass function at the low-mass 
end may provide a better agreement with observations for 
the low-z, intermediate mass galaxies, which can be tested 
by our future simulations with a larger volume. 

Although our new MVV wind model has new features 
and show a good agreement with observations, there are 
still many uncertainties in the details of the wind physics. 
For example, it is not clear when, where, and how the 
momentum-driven wind and the energy-driven wind domi- 
nates over one another. In the MVV wind models, we mixed 
the momentum- driven wind and the energy-driven wind 
based on the local gas density. This mixture method and 
criterion are based on a logical but an ad hoc argument, 
and they need to be carefully evaluated and extensively 
tested. From the comparisons with observations presented in 
this paper, it is still difficult to distinguish different models 
and identify a particularly favoured model. The overdensity- 
metaliicity relation (Figure [TT]) suggests that the EE run is 
the favoured model, but the Q(Civ) evolution shows that 
the EE run presents the worst agreement among different 
MVV runs. Despite of these uncertainties, the comparisons 
among the four MVV models provide some insight into the 
wind models. The result from EE run shows that it can 
reduce the SFR considerably, but this leads to significantly 
increased gas accretion at 2; < 4. It suggests that the energy- 
driven wind is a possibly legitimate mechanism for low mass 
and high- 2; galaxies, but a stronger outfiow model is required 
for high- mass, low- 2; galaxies. 

Also, the details of driving mechanism of the multiphase 
outfiow and its thermal composition are not well under- 
stood. For example, the observed values of t'wind and its 
functional dependence on various galaxy parameters, such 
as the stellar mass or SFR, are still unclear. We postulated 
'i^wind cx: (SFRy^^ based on Equations dH) and (|5]). How- 
eve r, it is well known that Vc oc mI^^ for early type galax- 
ies (Faber Jackson|ll976l ). where is the galaxy stellar 
mass. N oeske et al.l (|200^ showed a strong positive corre- 
lation of SFR oc M*. Combining the above two relation- 
ships suggests t'wind OC (SFR)^^^. In this paper, we adopt 
the relationship t;wind oc {SFRy^^ based on the arguments 
presented in Section [2^21 and also because some observations 
prefe r ^;wind oc (SFR)^^^ scaling (|Martinll2005l : IWeiner enil 
I2OO9I ). 

Given these current uncertainties in both theory and ob- 
servations, it is still somewhat premature to constrain and 
determine the best model among the four tested MVV wind 
models. Therefore, our new model may be yet another ap- 
proach to the phenomenologicai model of galactic wind in 
cosmologicai SPH simulations, but our goal is to develop 
a better model for star formation and feedback incremen- 
tally for galaxy formation simulations. To further discrim- 
inate the models, we intend to perform further analyses of 
our simulations with regard to the quasar absorption lines, 
e.g., Ly-a forest and damped Ly-a absorbers. At least the 
current comparisons between our simulations and observa- 
tions provide some improved understanding on the role of 
r] and Vwind in determining the cosmic SFR and IGM metal 
enrichment. 
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In summary, our new galactic outflow model, which re- 
flects the fact that the wind material can arise from different 
phases of ISM and captures the dependence of variable wind 
velocity on galaxy mass and SFR, provides better agree- 
ment with various observations, such as the cosmic SFR, 
gas overdensity-metallicity relationship, and Q(Civ). These 
comparisons suggest that the MVV wind model is definitely 
favoured over the CW model, and that it is a more viable 
choice as a wind model for the future simulations. 
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